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Abstract 



We establish the global existence and uniqueness of strong solutions to the initial 
boundary value problem for incompressible MHD equations in a bounded smooth 
domain of three spatial dimensions with initial density being allowed to have vac- 
uum, in particular, the initial density can vanish in a set of positive Lebessgue 
1 measure. More precisely, under the assumption that the production of the quan- 

tities || v / Po' Uo lll 2 (n) + 11-^0 |||a( n ) and ||V«o|li2 ( Q) + ||Vi7o|||2( n ) is suitably small, 
qq ! with the smallness depending only on the bound of the initial density and the do- 

] main, we prove that there is a unique strong solution to the Dirichlet problem of 

the incompressible MHD system. 
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Magnetohydrodynamics (MHD for short) is the study of the interaction between mag- 
netic fields and moving conducting fluids, which can be described by the following system 

p t + div(» = 0, (1.1) 

p(u t + (u ■ V)«) - pAu + Vp = (V x H) x H, (1.2) 

H t - XAH = V x (ux H), (1.3) 

divu = divH = 0, (1.4) 

where p = p(x,t) G M + denotes the density, u = u(x,t) G M 3 the fluid velocity, H = 
H(x, t) G M 3 the magnetic field, p = p(x, t) G R the pressure, positive constant p is called 
the kinematic viscosity and positive constant A is called the magnetic diffusivity. Usually, 
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we refer to equation (11. ip as the continuity equation, which represents the conservation 
lass of the mass, to ( II. 2p as the momentum conservation equation. Equation (11.31) is a 
considerably simplified version of the well known Maxiwell's equations (see [5j [lj |2]) by 
dropping Gauss's law and ignoring the displacement currents, and it's sometimes called the 
induction equation. As for the constraint divH = 0, it can be seen just as a restriction on 
the initial value H of H since divi^ = 0. Note that, using the condition ( II. 4p . equations 
(ll.2p and (jl.3p can be rewritten as 

p(u t + (u ■ V)u) - pAu + = (H ■ V)H, 
H t + (u- V)H - XAH = (H ■ V)u. 

There are a lot of literatures on the study of MHD. The global existence of weak 
solutions to the homogeneous incompressible MHD was proven in [3] and [I] long time 
ago, the density dependent case was later proven in [5] by using Lions's method [6]. The 
corresponding results on the global existence of weak solutions to the compressible MHD 
are proven in [HI [9j [131 HQ] by using the method exploited in [7] (see also [HI [12]), where in 
the first paper the isentropic case is considered, while in the other three ones concern the 
non-isothermal model. Local existence and uniqueness of strong solutions can be found 
in P3] and [15] for incompressible model, in [T7] and [18] for the compressible model, 
where in the last three papers the non-isothermal model are considered. Global existence 
and uniqueness of strong solutions to the incompressible MHD with vacuum in two space 
dimensions is proved in [16], recently, and this work is new even for the inhomogeneous 
incompressible Navier-Stokes equations. While the three dimensional case, with small 
initial data and away from vacuum, is proven in |19j . 

One of the most important questions for system (jl.lj) — (jl.4p is to prove the global 
existence and uniqueness of solutions (p,u,p) satisfying the initial condition 

(p,u,H)\ t =o = (po,u ,H ) in fi, (1.5) 

and boundary condition 

u\ dQ = 0, H-n\ d Q = 0, (V x H) x n\ dn = 0, (1.6) 

where n is the unit outward norm vector on dQ. 

The aim of this paper is to prove the global existence and uniqueness of strong solutions 
to system (ll.ll) - (ll.6p with the initial data being allowed to have vacuum. For 1 < p < oo, 
we denote by H^Hp the L P (Q) norm of the function u. The definition of strong solution is 
stated in the following: 

Definition 1.1. We call (p, u, H) a strong solution to the system U.l\) - I[l~b}) on (0, T) 

if(p, u, H) satisfies U.1\) - [L4\ ) a.e. in fix (0, T) for some pressure function p, and satisfies 
the initial condition U.5\) and boundary condition (HE) , with the regularity 

p E L°°(Q T ) n L°°(0, T; H\n)), Pt E L°°(0, T; L 2 (fi)), 

u e l°°(o, T; H*(n) n H 2 (tt)) n l 2 (o, T-, w 2fi (n)), 

H G L°°(0, T; H 2 (n)) n L 2 (0, T; W 2 ' 6 (n)), 

u t ,H t G L 2 (0,T; H\n)), ^pu t ,H t G L°°(0, T; L 2 (fi)). 



Our main result is stated bellow: 

Theorem 1.1. Let Q be a bounded smooth domain in M 3 and p a positive number. 
Assume that the initial data [pq,Uq) satisfies the conditions 

< po < p, po e H\n), u e Hl(tt) n H 2 (n), H e H 2 (tt) 

and the compatibility condition 

divu = divH = 0, H ■ n\ dn = 0, (V x H ) x n\ dn = 0, 
Au + (H -V)H - Vp = VPo9o in 

for some {po,9o) £ H l (Q) x L 2 (Q). Then there exists a positive constant So depending 
only on p and Q, such that, if 

(\\Vp~ou \\1 + \\H \\l)(\\Vu \\l + llVHolU) < s 0) (1.7) 

then the system ( fl.ij) -( TOl) has a unique global strong solution. 

Remark 1.1. The quantity (H-v/p^lli + II-^IIDGI^HII + 11^-^111) ^ s sca ^ n 9 invariable 
under the scaling transform 

px(x,t) =p(\x,X 2 t), p x (x,t) = \ 2 p(\x, \ 2 t), 
u\(x, t) = \u(\x, X 2 t), H\(x, t) = \H(\x, X 2 t), 

and thus Theorem \l.l\ can be viewed as a result on the global existence of strong solutions 
with vacuum in critical space. 



2 Proof of Theorem 11.1 

Throughout this section, we denote 

1 1 -x/Po^o 1 1 1 + ll-^o H2 = Co- 
Definition 2.1. A finite time T* is called the finite blow-up time if 
$(T) < 00 for all < T < T* and lim $(T) = oo, 

where the function $(T) is given by 

$(T) = sup (||(Vp,a)|| 2 + \\(u,H)\\ H 2 + \\H t \\ 2 + \\^u t \\ 2 ) 

0<t<T 
cT 



+ [ (\\(u,H)\\ 2 w2fi + \\(u tl H t )fm)dt. 
Jo 



We will use the following lemma, which states the local existence and blow up criterion 
of local strong solutions. 



Lemma 2.1. (See [14]) Under the conditions of Theorem ! 1. 11 there hold 

(i) (Local existence) there exists a small time T* and a unique strong solution on 
(0,T*) ; 

(ii) (Blow-up criterion) T* is the finite blow-up time of (p,u, H) if and only if 

[ \\(Vu, VH)\\ 8 2 dt < oo for all < T < T„ and [ \\{Vu, VH)\\ 8 2 dt = oo. 
Jo Jo 

To proof the global existence of strong solutions, we need to extend the local strong 
solution given in the above lemma to be a global one. For this purpose, we need do some 
a priori estimates on the local strong solutions. The following two lemmas give the energy 
estimates on the strong solutions, where the first one concerns the basic energy estimates, 
while the second one concerns the higher order estimates. 

Lemma 2.2. Let (p,u,H) be a strong solution to system ( fi.-/j) -( TO)) on (0, T). Then, 
there holds 

*t 



(\\^pu\\l+\\H\\l)(t)+2 / (||V«|| 2 +\\VH\\l)ds=\\vfau \\l+\\H { 

Jo 



o\\l 



for any < t < T. 

Proof. Multiply (11. 2\\ by u and integrate over Q, by the aid of (11.11) . we obtain after 
integration by parts that 

f ?-\u\ 2 dx + [ \Vu\ 2 dx= [ (H ■ VH) ■ udx = — f (H ■ V)w • Hdx. (2.1) 
dt Jn 2 J n J n J n 

Multiply (ll.3p by H and integrate over Q, we obtain after integration by parts that 

4 I ^—dx+ [ \VH\ 2 dx= [ (H-V)u-Hdx. (2.2) 
dt Jn 2 Jn Jn 

Summing ( 12. ip with (I2.2p up, and integrating the resulting equation with respect to t, we 
obtain 

f (p\u\ 2 + \H\ 2 )dx + 2 f [ (|Vw| 2 + \VH\ 2 )dxds = [ (p \u \ 2 + \H \ 2 )dx, 
Jn Jo Jn Jn 

completing the proof. □ 

Lemma 2.3. Let (p,u,H) be a strong solution to system $l.l\) - [T7b}) on (0,T). Then, 
there holds 

t [\\V 2 u\\ 2 2 + \\V 2 H\\ 2 2 )ds 



sup(\\Vu\\ 2 2 + \\VH\\j)+ [ 

0<s<t Jo 

<2(||Vmo||1 + ||VJy ||i) + C sup (C 2 ||Vm||^ + C \\VH\\ 3 2 ) 

0<s<t 

+ CC sup (||Vtz|| 2 + ||Vif || 2 ) / (||V 2 m||^ + ||V 2 if||^)^ 

0<s<t Jo 

for any t G (0, T). 



Proof. Multiplying (11. ip by Ut and integration by parts yields 

d f |Vu' 2 
dt 



dx + p\u t \ 2 dx 
Jn 

/ [(H ■ V)H ■ ut - p{u ■ V)u ■ u t ]dx = - [p(u ■ V)u ■ u t + (H ■ V)u t ■ H]dx 
Jn Jn 

- f(H- V)« ■ Hdx + [ [{H t -V)u-H+{H- V)u ■ H t - p{u ■ V)u ■ u t ]dx 
dt Jn Jn 

{H-V)u-Hdx + j {^p\u t \ 2 + h [ H t \ 2 ^j dx + A j \\H\ 2 \Vu\ 2 + p\u\ 2 \Vu\ 2 )dx, 



d 

< 

~ dt 



and thus 

d f f\Vu\ 



+ (H ■ V)u ■ Hj dx + j p\u t \ 2 dx 



< f (\p\ut\ 2 + \\H t \ 2 ) dx + A I (\H\ 2 \Vu\ 2 + p\u\ 2 \Vu\ 2 )dx. (2.3) 



.2' 1 1 4 

Multiplying (II. 3p by H t and integration by parts yields 



d f \V^ dx+ f \ Ht \2 dx= f {[H-V)u-{u-V)H]H t dx 
dt Jq 2 Jq Jq 

<\ j \H t \ 2 dx+ [ (\H\ 2 \Vu\ 2 + \u\ 2 \VH\ 2 )dx. (2.4) 



"4.... 

Summing (12.31) with fl2.4jl up, it follows 

d 
dt 



and thus 



(|V«r + \VH\ 2 + 2(H ■ V)u ■ H)dx + / (p\u t \ 2 + \H t \ 2 )dx 
i Jn 

<10 f (\H\ 2 \Vu\ 2 + \u\ 2 \VH\ 2 + p\u\ 2 \Vu\ 2 )dx, 
Jn 

sup (\\Vu\\l + \\VH\\l) + [\\\^u t \\ 2 2 + \\H t \\l)ds 

)<s<t Jq 



0<s<t 

<(||Vu ||2+ l|V^o||2)+4 sup / \H\ 2 \Vu\dx 

o<s<tJn 

+ 10 f [ {\H\ 2 \Vu\ 2 + \u\ 2 \VH\ 2 + p\u\ 2 \Vu\ 2 )dxds. (2.5) 
Jo Jn 

Applying H 2 estimates to Stokes equations and elliptic equations, it follows from (II. 2p 
and (I1.3P that 

\\V 2 u\\l <C(\\pu t \\ 2 2 + \\p(u ■ V)u\\l + \\(H- V)H\\l) 
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<C{\\Jpu t \\l + \\^p(u • V)«||! + \\(H- V)H\\l) (2.6) 

and 

|| V 2 ^||l < C(\\H t \\l + \\(u. V)H\\l + \\(H. V)u\\l), (2.7) 

where we have used the fact that < p < p. Combining (12.51) with (12. 6p . together with 
(12. 7p . we obtain 

sup(||v«||! + ||wr||2)+ / (ll^ll' + liv^ + lliytll^ + nv 2 ^!^)^ 

0<s<t Jo 

<(l|V«o||2+ || Vi?o||l) +4 sup I \H\ 2 \Vu\dx 

o<s<tJn 

+ C f f [|#| 2 (|Vu| 2 + \VH\ 2 ) + \u\ 2 \VH\ 2 + p\u\ 2 \Vu\ 2 }dxds. (2.8) 
Jo Jn 

It follows from Holder inequality, Sobolev inequality and Young inequality that 

\H\ 2 \Vu\dx <||#|| 2 ||V M || 2 < \\H \\l /2 \\H\\T \\Vu\\ 2 

<C||#|| 2 /2 ||Vif ||5 /2 ||V«|| a < s\\Vu\\ 2 2 + CII^HallVfTHl, 

P \u\ 2 \vu\ 2 dx <c||vHWMl6l|Vu|||j < cii^iisiiv^ihiiv^ii 2 , 

n 

\u\ 2 \VH\ 2 dx <||ii|| 2 ||V J f/|| 2 < C||V W || 2 ||Vtf || 2 ||Vtf || 6 

<C||V M || 2 ||VH|| 2 ||V 2 iy|| 2 < e\\V 2 H\\ 2 2 + CHV^HVtf || 2 , 



J \H\ 2 (\Vu\ 2 + \VH\ 2 )de <||if|| 2 ||#|| 6 (|| Vti|| 2 + \\vh\\D 

<C\\H\\ 2 \\VHU\\V 2 u\\ 2 2 + \\V 2 H\\ 2 2 ). 

Substituting the above inequalities into (I2.8p . taking e small enough, it follows from 
Lemma 12.21 that 

sup (||V W || 2 + || Vtf|| 2 ) + AhVp^II! + II + ll^tlla + HV^lDda 

0<s<t Jo 



<2(||V« || 2 + HVffolli) + C sup ||iJ|| 2 ||V#||i + C sup (||V^|| 2 ||Vw|| 2 + \\H\\ 2 \\VH\\ 2 ) 

0<s<t 0<s<t 

x A||V 2 m|| 2 + ||V 2 //|| 2 )c/s + C sup ||Vu||2 f \\VH\\lds 

JO 0<s<t Jo 

<2(\\Vu \\t + \\VHo\\l) + C sup (C^IIV^HI + C ||Vi?||l) 

0<s<* 

+ CC sup(||Vw|| 2 + ||Vif|| 2 ) / (||V 2 W || 2 + ||V 2 #|| 2 )cfe, 

0<s<t Jo 

completing the proof. □ 
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By the aid of the above two lemmas, we can prove the following a priori estimates. 

Lemma 2.4. Let (p,u,H) be a strong solution to system ( fi. ij) -( TO)) on (0, T). Then, 
there exists a positive constant e depending only on p and Q, such that 

snp (||Vu||jj + + [ (\\V 2 u\\l + \\^ H \\l)dt < K\\Vu Q \\l + HVflollI), 

0<t<T Jo 

provided 

C (\\Vu \\l+\\VH \\l)<e . 
Proof. Define functions E(t) and <£>(£) as follows 

E(t)= sup(\\Vu\\l+\\VH\\ 2 2 )+ [ (||V 2 u||l + ||V 2 #||l)ds, 

0<s<t Jo 

m = Cl sup(||Vn|| 2 + ||V^|| 2 ). 

0<s<t 

In view of the regularities of u and H, one can easily check that both E(t) and $(£) are 
continuous functions on [0,T]. By Lemma [2.31 there is a positive constant C*, such that 

E(t) < 2(\\Vu § + \\VHoWi) + C^mf' 2 + mW)- (2-9) 

We take 

Co = min 



and suppose that 
We claim that 



1 1 

32a' 128C 2 f 

^o(l|V«o||2+l|V^o||2)<^0. 



^Xmm^,^^, 0<t<T. 



Otherwise, by the continuity and monotonicity of $(£), there is T G (0,T], such that 

$(T ) = min(-^,-l-V (2.10) 



.4C*'16C 2 , 
On account of OTTO]) , it follows from (P23D that 

£(T ) < 2(||V^o||l + IIVfTolll) + \e(T ), 

and thus 

S(ro)<4(||Vuo||^ + ||VH ||i). 
Recalling the definition of E(t) and $(t), we deduce from the above inequality that 



$(T ) <C££(T ) < 4C 2 (||Vu ||S + HV^olla) < = min 



1 1 



8C* ' 32C 2 j ' 
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contradicting to (12.101) . This contradiction implies that the claim is true. 
By the aid of the claim we proved in the above, it follows from (I2.9P that 

E{t)<4(\\Vuo\\l + \\VH \\l), 0<t<T, 

completing the proof. □ 

Now, we can give the proof of Theorem 11.11 

Proof of Theorem ll.lt Let £ be the constant stated in Lemma 12.41 and suppose 
that the initial data satisfies 

(llV^ollI + ||i?oll2)(l|Vtio||! + HVfloll!) < e . 

By Lemma [27T] there is a unique strong solution (p,u,H) to system (ll.ll) - (ll.6l) . Extend 
such local solution to the maximal existence time interval [0,T*). We will prove that 
T* = oo. Suppose, by contradiction, that < oo. By Lemma 12.11 the time T* can be 
characterized as follows 

/ ||(Vu,V#)|||<ft<oo, for all < T < T„ 
Jo 

and 

/ * \\(Vu,VH)\\*dt = oo. (2.11) 
Jo 

By Lemma [2.4[ f° r an Y < T < T*, there holds 

sup (\\Vu\\l + ||ViT||=) < ±(\\Vu f 2 + \\VH \\l), 

0<t<T 

and thus 

sup (\\Vu\\l + \\VH\\l) < A(\\Vu \\ 2 2 + IIVHolU), 

0<t<T» 

which implies 

/ * \\{Vu,VH)\\ 8 2 dt < 4(||V«o||l + l|VFo||^)T* < oo, 

contradicting to (12.111) . This contradiction provides us that T* = oo, and thus we obtain 
a global strong solution. The proof is complete. 
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